ON SOME EQUATIONS THAT CAN BE BROUGHT 
TO THE EQUATIONS OF DIFFUSION TYPE. 
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Abstract. For the system of second order quasilinear parabolic equations the problem 
of reducing them to the equations of diffusion type is considered. In non-degenerate 
case an effective algorithm for solving this problem is suggested. 



1. Introduction. 

Let's consider the system of differential equations defined by some matrix A = 
A^y 1 , . . . ,y n ) with non-zero determinant: 

3=1 

The equations (1.1) do not form the set invariant under the point transformations 
given by the following changes of variables: 

(1-2) 

y n = y 1 (y\... ,y n ). 

However, one can expand this set up to the point-invariant class of equations: 
dy^ 1 _^ f Q 1 ^ 1 ^ ^ Qy^ Qy^ \ 

(L3) Tt = ^ A H 9? + EE r «& -^ J • 

j=l \ r=l a=l / 

Here coefficients A* and T 3 rs depend on y 1 , ... , y n so that det A ^ 0. Denote by 
T the Jacoby matrix for the change of variables (1.2) and denote by S the Jacoby 
matrix for inverse change of variables: 

g i = djf_ T , = dff_ 

3 dy 3 ' 3 dy 3 ' 

In geometry matrices S and T are called direct and inverse transition matrices 
respectively (this implies that transition from y , . . . , y n to y , . . . , y n is consid- 
ered as direct transition). 
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Theorem 1.1. Under the point transformation (1.2) the equations (1.3) transfer to 
the equations of the same form 

HA) ^ = fiite + ffw djTdyA 

1 ■ ' dt ^ 3 \ dx 2 rs dx dx )' 

] = l \ r=l s=l / 

Parameters A % -, T J rs , A % - and r 3 rs in (1.3) and (1.4) are related as follows 

n n 

(1.5) A*=Y,Y, S rn T ?A™, 

rn—l p—1 

TL 7X TL TL C\ r f^'(Tl 

(i-6) = E E E s ™ T f T l f m + E s " 1 ' 



m—lp—lq—l m—1 



If wc suppose y 1 , . . . , y n to be local coordinates on some manifold M, then due to 
theorem 1.1 parameters A\ define a tensor field A of the type (1, 1) on M . Parameters 
r^. in turn can be interpreted as components of some affine connection T on M. This 
enables us to apply various geometrical methods to the investigation of the equations 
(1.3). Here we consider only one problem due to equations (1.3), which appears to 
have graceful geometrical solution. Similar problems for other classes of differential 
equations were considered in [1-6]. 

2. Reduction to the equations of diffusion type. 

Equations (1.1) arise in various models describing the diffusion phenomena in mul- 
ticomponent mixtures (see [7-24]) and in theory of integrability (see [25-30]). They 
are called the equations of diffusion type. As for the equations (1.3) those of them, 
that can be brought to the form (1.1) by means of point transformation (1.2), should 
obviously be involved to the class of equations of diffusion type. This gives rise to the 
problem of finding effective description for this class of equations. 

Problem. What condition for tensor field A and affine connection T should be ful- 
hllcd on M in order to guarantee the existence of point transformation (1.2) that 
brings equations (1.2) to the diffusion form (1.1)? 

Suppose that one could bring (1.3) to the form (1.1) in local variables y 1 , ... , y n . 
Then in such variables for components of affine connection T we have 




Here B™ are the components of the matrix B which is inverse to the matrix A formed 
by components of operator field A in local coordinates y 1 , ... , y n . Now let's apply 
formulas (1.5) in order to transform the components of the field A in (2.1): 



n n 

A s — V* \^ T s <>* A r 
r—1 i—1 
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Further differentiate this relationship with respect to y q : 

(") § = 1 1 (f * « + 1 § « 

Let's express the derivatives dA r i /dy q , dS p /dy q , and dT^/dy q through partial deriva- 
tives in initial variables y 1 , ... , y n . We shall do this by means of formulas 

OAl = ^ qj dAl dJl = ^ jdTl d_Sl » » qj dTl qa 

Qyq ^1 Qyj ' Qyq °9 Qyj ' 0» 2^ ^ ^ ^ ^ ■ 

Substitute them into (2.2), then substitute the resulting expression into the formula 
(2.1) for the components of T. Missing expression for the derivative dA q /dy p we get 
by transposing indices p and q: 

Tt Tt TL Tt i-lT^S 

r — 1 s— 1 i—l j — 1 

n n n n 



1 '\ '\ '\ '\ _ 9^ 

9 X! X! ^2 ^2 T r -faj i S l S q + S% q S p) - 

r—1 s—1 i — 1 j—1 



2^^^^ s r dyi 

s=l j=l a=l {3=1 



1 n n n n nmrt 



Now substitute the above expression for r™ into the formula (1.6). As a result we 
obtain the following formula for Y\- : 

n R*: / F> Ar f)A r \ n n n tp <rq , ap cq 

r=l v y y 7 r =lp=lg=l 

Here 5? is Kronecker's delta-symbol (representing components of unit matrix): 

1 1 for q ^ i. 

By we denote the components of operator field B = A -1 . The entrance of matrices 
S and T into (2.3) is completely defined by the quantities 9 r pq : 

(9 d) ft r = <? r — E = f r 

y ■ ' PI £^ J ™ Q q 2^ °™ Q p Q g ■ 

m=l a m=l a a 

From (2.4) we see that 9 pq = 9 qp . If we know that system of equations (1.3) admits the 
coordinates y 1 , . . . , y n , where it takes diffusion form (1.1), and if these coordinates 
are given, then by means of formula (2.4) we derive 9 pq for further substitution into 
the relationship (2.3). 
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Otherwise, when coordinates y , . . . , y n are not given, one should determine Q r v 



from (2.3) and should consider (2.4) as the equations that determine required change 
of variables (1.2). We write these equations as 

Q r J lm n 
y r=l 

Equations (2.5) form complete system of linear Pfaff equations respective to the com- 
ponents of the matrix T. They are well known in differential geometry. Compatibility 
condition for (2.5) is derived from permutability of partial derivatives: 

, . d 2 T™ Q 2 Tl n 
(2.6) S— 



Qyq Qyk Qyk Qyq 

Let's calculate both sides of (2.6) according to equations (2.5). As a result we get the 
relationship, which is known as the condition of "zero curvature" 1 : 



d9Z 89" 

c + 

r— 1 r— 1 



fr, rj\ kp qp \ " nm nr \ ~* nm nr 

y 2J ) -QyT - ~QyT + 1^ °qr »kp - kr V qp - «■ 



On the base of theory of Pfaff equations (see, for instance, appendix A in [4]) one can 
prove the following theorem. 

Theorem 2.1. Let v be some point of the manifold M. Matrix Pfaff equations (2.5) 
have solution with non-zero determinant det T ^ in some neighborhood of the point 
v on M if and only if the condition (2.7) is fulhlled. 

Quantities Q r vq are symmetric in lower pair of indices p and q. Hence for arbitrary 
solution of the equations (2.5) we have 

<9T D m QT™ 
(2.8) ; - 



Q y q QyP ' 

This condition (2.8) in turn is the compatibility condition for the equations 

Qy" 



dyP 



T™(y\...,y n ) 



which determine new variables y 1 , ... , y n . We shall formulate the result of the above 
speculations in form of two theorems. 

Theorem 2.2. Quantities d r pq = r qp are defined by some change of local coordinates 
in some neighborhood of the point v on M according to the formulas (2.4) if and only 
if the conditions of "zero curvature" (2.7) are fulhlled. 



^^Note that "zero curvature" condition (2.7) docs not mean that the curvature tensor of the afBnc 
connection T is zero. 
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Theorem 2.3. System of equations (1.3) can be brought to the diffusion form (1.1) 
by some change of local coordinates in some neighborhood of the point v on M if and 
only if one can find quantities d r vq = 9 qp satisfying both conditions (2.3) and (2.7) 
simultaneously. 

Theorem 2.3 distinguishes two cases in the theory of quasilinear equations of the 
form (1.3): first is the case of general position and second is the degenerate 
case. In the case of general position quantities Q r pq are uniquely determined by the 
relationships (2.3). Substituting them into (2.7) further we can get required equations 
for A and T which guarantee that the equations (1.3) are reducible to the form (1.1). 

In degenerate case relationships (2.3) do not determine the quantities Q r pq at all 
or determine with some extent of uncertainty (not uniquely). In this case we should 
analyze the relationships (2.3) for their solvability with respect to Q r pq and for their 
compatibility with (2.7). 

3. Non-degeneracy condition. 

The non-degeneracy condition distinguishing the case of general position, as formu- 
lated above, is the condition of unique solvability of the equations (2.3) with respect 
to Q r . These equations can be written as follows: 

A^S q 4 +A P 5j A .„_, i (dA- OA- 



p=lq=l k=l V W 



The equations (3.1) are linear equations respective to d r pq . Note that they breaks into 
n separate parts. For each r = 1, ... , n we have closed system of n 2 linear equations. 
Matrices of all these systems are the same, they are defined by the components of 
tensor field A in the following form: 

(3 . 2) ^ = ^q + Atq + ASq + AW 

Denote by w\- right hand sides of (3.1). Then these equations can be written as 

n n 

(3-3) EE A «^ = <- 

p=l q=l 

Let v be some fixed point on the manifold M . Denote by V the tangent space to 
M at this point: V = T V (M). By V* denote the dual space for V. Then let's consider 
a tensor product of two samples of dual space: 

W = V* ® V* = T$(v,M). 

According to this definition W is a space of bilinear forms on the space V. It is 
represented as a direct sum of two subspaces: 

(3.4) W = W sym © W skow . 



First is the subspace of symmetric bilinear forms and second is the subspace of skew- 
symmetric forms. Tensor A can be treated as a linear operator in the space W. Both 
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subspaces in the sum (3.4) are invariant with respect to the action of the operator 
A: W — ► W, restriction of A to M^kew being identically zero (since A^ 9 is symmetric 
in both pairs of indices upper and lower). 

Right hand side of (3.3) is symmetric in i and j, quantities Q r that should be 
determined from the equations (3.1), are also symmetric in their lower indices p and 
q. Therefore the above condition determining the case of general position can be 
formulated as the condition of non-degeneracy for the restriction of the operator A 
to the subspace of symmetric bilinear forms: 

^sym • l^sym * l^sym ■ 

In other words this is written as 

(3.5) detA sym ^0. 

Components of the operator A are uniquely determined by the components of A 
according to the formula (3.2). But nevertheless (3.5) doesn't follow immediately 
from det A ^ 0. One should check it up separately. 

4. Effectivization of the non-degeneracy condition. 

Let's study the question of how to check up the condition (3.5), if the matrix 
of operator field A is given. We shall start with the case when operator A has 
purely simple spectrum, i. e. when its eigenvalues Ai, ... , A„ are distinct and its 
characteristic polynomial has no multiple roots: 

n n 

(4.1) /(A) = dct(A - AI) = (-A)" + ^> s (n) (-A)"- s = J](A S - A). 

8=1 S=l 

Here g\ (n), ... , a n (n) are basic symmetric polynomials defined by formula 

(4.2) a s (n)= ]T...]T A, • . . . • A, . 

l^«i < .. . <i s 

Last of these polynomials is equal to the determinant of the operator A (see more 
details in [31]). 

Operator A with simple spectrum in finite dimensional space V = T V {M) is di- 
agonalizable, it has a base of eigenvectors ei, ... , e„. However, the construction of 
such base may require complexification of the space V, since roots or characteristic 
polynomial (4.1) are complex numbers in general. Denote by h 1 , ... , h™ the dual 
base in dual space V*. We have the following relationships: 

(4.3) Ae, = A 4 e 4 , A*h* = A,h l , h l (e k ) = 5\. 

Here A* : V* — > V* is a conjugate operator for A. By means of covectors h 1 , . . . , h" 
we shall build the base in the space of symmetric bilinear forms W syra : 

(4.4) w ij = h* ® h j + h J ® h\ where 1 < i ^ j < n. 
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Denote by N the dimension of the space W syra . On counting the number of vectors 
in the base (4.4) we find 

(4.5) N = dimWsyni = !^±ll. 

Due to (3.2) and (4.3) the result of applying A to the vector h p ® h 9 is given by 

(4.6) A(h p ®h q ) = A P + A| ? h p ®h q . 

Now let's apply the operator (4.6) to the vectors of the base (4.4). By means of direct 
calculations we get 

(4.7) A sym (w«) = A(w«) = w«. 

The relationship (4.7) means that operator A sym is diagonalizable, base (4.4) be- 
ing the base of eigenvectors of this operator. Eigenvalues mj corresponding to the 
eigenvectors w y are given by the formula 

(4.8) Mtj = A ' + 2 — , where 1 i sC j < n. 
Using (4.8) we can calculate the determinant of the operator A S y m : 

(4.9) det A sym = a N (N) = J] (hl±hl 

Let's introduce lexicographic ordering in the set of pairs of indices where 
i ^ j. Say that (ii,ji) < («2, J2) if ji < J2 or if «i < i 2 in the case of coincidence 
ji = ]2 ■ Then consider the sum 

(4.10) a s ( N )= E-E (Atl+A3l), 2 ;' (Als+AjJ - 

(ii,ji)< ... <(i s ,js) 

Formula (4.10) determines basic symmetric polynomials for the set of numbers fiij 
from (4.8). For s — N it is equivalent to (4.9). Its easy to note, that right hand sides 
of (4.10) do not change under any transposition of numbers Ai, ... , A„. They are 
symmetric polynomials in Ai, ... , A„. Hence according the well known theorem (see 
[31], chapter 6, §2) one can express them trough basic symmetric polynomials: 

(4.11) a s (N) = P s (a 1 (n),... ,a n (n)), where s = l,...,N. 

For any particular n we can derive explicit expressions for polynomials P s (a\, ... ,<r n ). 
For instance, if n = 2, we have 



(4.12) Pl = l a ^ P2 = \o\ + °2, Pz = \w\, 
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From (4.5) here we get N = 3. Since 0-2(2) = det^4 is non-zero, for n = 2 the 
non-degeneracy condition (3.5) can be written as 

(4.13) o-i (2) = tr A ^ 0. 

Theorem 4.1. For n = 2 the system of equations (1.3) belongs to the non-degenerate 
class if and only if det A^O and tr A ^ 0. 

The relationship (4.13) means that operator A shouldn't have two eigenvalues 
which differ only in sign. In this form the non-degeneracy condition can be stated in 
multidimensional case as well. For arbitrary n it follows from (4.9). Indeed, we have 

(4-i4) n(^)=de tA .n(^Wo, 

which is equivalent to det A ^ and \ ^ —Xj for i ^ j. 

In order to check up the condition (4.14) one need not find eigenvalues of the 
operator A. This condition can be written in terms of coefficients of characteristic 
polynomial (4.1), provided we have explicit expressions for polynomials in (4.11). The 
latter is algorithmically solvable problem. 

Another effective way of checking up the condition (4.14) is based on the theory 
of resultants. Let /(A) be a characteristic polynomial for the operator A. Its roots 
are non-zero since det A 7^ 0, therefore the condition A, ^ — Xj for i ^ j means that 
polynomials /(A) and /(—A) have no common roots. Then resultant of these two 
polynomials is non-zero (see [31] or [32]): 

(4.15) Res[det(A- AI),det(A + A I)] ^0. 

Theorem 4.2. System of equations (1.3) belongs to the non-degenerate class if and 
only if det A^O and condition (4.15) is fulhlled. 

5. Finding inverse operator for A sym . 

Suppose that non-degeneracy condition stated in theorem 4.2 is fulfilled. Let's find 
regular way for inverting the operator A sym : W sym — > W sym . One need to do it for 
to find B T pq from the equations (3.3). Let /(A) be a characteristic polynomial for the 
operator A. Lets consider pair of polynomials /(A + fi) and f(fi — A), where fi is a 
parameter. These polynomials are factored as follows: 

n 

/(A + M ) = (-irJ](A-(A i - M )) 

(5-1) „ 4=1 

/(/i-A) = IJ(A-(/i-A i )). 

i=l 

For to calculate the resultant of polynomials (5.1) let's apply the following theorem 
of algebra (see [31] or [32] for proof). 
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Theorem 5.1. Suppose that two polynomials g(X) and h(X) with leading coefficients 
a and b are factored to a product of linear terms 

g(X) = a-{X-a 1 )-...-(X-a n ), 
{ '' h(X) = b-(X-/3 1 )-...-(X-/3 m ). 

Then resultant of these polynomials is given by the formula 

n m 

(5.3) Res [ 5 (A), MA)] = a" b m JJ JJ( tti - ft). 

i=1.7'=l 



Comparing (5.1) with (5.2) from (5.3) we obtain the resultant of polynomials (5.1): 

n n 

Resf/Ou - A), /(A + fj,)] = (-1)" 1] Il( 2 ^ - A < - A ^') = 

i=l j=l 

(5.4) 

= 2 " 2 n(^-M)n(^-^ 

Let's compare (5.4) with formula (4.8) for eigenvalues of the operator A S y m . This 
comparison gives the formula 

(5.5) ^( M ) 2 = 2-" 2 f[p) Rcs[/(/i - A), /(A + rf\ . 

Here ip(p) = det(A sym — /ul) is a characteristic polynomial for the operator A sym . 
Its remarkable that formula (5.5) express the square of characteristic polynomial of 
operator A sym in terms of characteristic polynomial of operator A. 

Denote by e(fi) the polynomial in right hand side of (5.5). According to the well 
known theorem of Hamilton and Cayley (see [32] or [33] ) when substituting the oper- 
ator A sym into its characteristic polynomial (p(fi) we should get zero operator. This 
is true for the polynomial e(p) = </?(/x) 2 as well: 

M 

(5.6) e(A sym ) = ^ £j A: ym = 0. 

i=0 

Here M — 2 N = n(n + 1) is the degree of polynomial e(p) = (fi(^) 2 ■ Now we can 
rewrite (5.6) as follows 

(5-7) (j^* A ^") ' Asym = ~ £ ° 1 

Quantity e in right hand side of (5.7) is the value of polynomial e(A) at A = 0: 

so = e(0) = <^(0) 2 = (detA sym ) 2 . 
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Due to non-degeneracy condition (3.5) the quantity £o is non-zero. Therefore we can 
use (5.7) in order to calculate the operator A^^: 



M 



A- 1 =-V 

sym / , 



i e ° 



We can expand operator : W sym — > W sym to be identically zero in the subspace of 
skew-symmetric bilinear forms Wskcw Though being degenerate, expanded operator 
A~ y \ n is defined in the whole space of bilinear forms W. One can define this operator 
by the following formula 

M 

(5.8) Ai = -Er rl - 

Right hand side of (5.8) is interpreted as a tensor field D of the type (2, 2) whose 
components D* J are symmetric in upper and lower indices. 

Theorem 5.2. If detA ^ and if non-degeneracy condition (3.5) written in the 
form (4.15) is fulfilled, then equations (3.3) can be solved with respect to r vq : 



(5-9) ^ = EE D 

1=1 ] = 1 



ii r 
W- ■ 
pq ^i] 



Here D l J q are components of tensor field D = A^^ defined by the formula (5.8). 

Naturally we can choose more straightforward way of calculating components of 
tensor field D other than formula (5.8). In order to do it we should only remember 
that D defines an operator inverse to the operator A sym in subspace W sym : 

(5.10) D • A sym = I. 

Written in coordinates the relationship (5.10) looks like 

i=l j=l 2 

This is the system of N 2 linear equations with respect to D % J S , where ./V is defined by 
(4.5). Non-degeneracy condition (4.15) guarantees the compatibility of these equa- 
tions and the uniqueness of their common solution. 

6. Criterion for being in diffusion class. 

On solving the equations (3.3) respective to Q r in the form (5.9) we can substitute 
this expression into (2.7) and we can apply theorem 2.3. 
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Theorem 6.1. In non-degenerate case equations (1.3) can be brought to the diffusion 
form (1.1) by means of some point transformation (1.2) if and only if the quantities 

do satisfy "zero curvature" condition (2.7). 

Due to formulas derived in sections 3, 4 and 5 we can check up non-degeneracy 
condition and we can explicitly calculate D l J s in those variables y 1 , ... , y n , where 
the components of operator field A are initially given. Therefore theorem 2.3 is 
an effective criterion for testing if the equation belong to the diffusion class in non- 
degenerate case. The algorithm given by this theorem can be easily realized by means 
of any of now existing program packages for symbolic calculations. 

Note. All results in sections 3, 4 and 5 were obtained under the assumption that 
operator A has purely simple spectrum. However, they remain true for arbitrary 
operators, since any operator with degenerate spectrum can be obtained as a limit of 
some sequence of operators with purely simple spectrum. 

7. Explicit formulas for n = 2. 

As it was shown above, in two-dimensional case the non-degeneracy condition 
reduces to the following relationships 



detA^O, trA^O 

(see theorem 4.1). Using (4.12) we can find explicit expression for the characteristic 
polynomial ip(u) of the operator A sym . Its coefficients are expressed in terms of 
parameters o\ = tr A and <r 2 = dot A: 

(7-1) ^( M ) = - M 3 + — M 2 -^ ^+^-- 

From operator equality ip(A syr[L ) = due to (7.1) we get the following expression: 
(7.2) A- = 2A 2 -3. 1 A + (^ + 2, 2 )I^ 

CTl 0"2 

Substituting (3.2) into (7.2) we obtain the components of the tensor field D = A~^: 

n ll _ ^2 + {AD 2 22 _ 2 + (A\) 2 

01 CT 2 <J\ (72 

n ll _ (A-l) 2 n 22 _ (^l) 2 

22 ~~ „ „ ' ~~ „ „ ' 

01 02 01 02 

4 2 4 2 A 1 A 1 

n 12 _ n 21 _ _£M_£*2 n 12 _ n 21 _ A 2 71 1 

u \\ — u \\ — ; -^22 ~ ^22 — ! 

01 02 01 02 

A 1 A 2 A 2 A 1 

n 11 — n 11 — 2 2 n 22 — n 22 — 1 1 

^12 — u 2i — > u \2 — u 2\ — ' 

01 02 0i 02 
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A 1 A 2 A 1 A 2 

n l2 _ n 2l _ A l ^2 n l2 _ n 2l _ ^2 

u 12 ~ U Y1 ~ J u 1\ ~ U 1\ — • 

(T\(T2 <J\Oi 



Now we should substitute the above expressions into (6.1) for to find parameters 
9p q and substitute them into (2.7). As a result we obtain four differential equations 
binding Aj and Y^. If they hold, then the equations (1.3) can be brought to the form 
(1.1) in two-dimensional case n = 2. 
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